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1 Introduction 

This paper is concerned with the existence of maximal, asymptotically flat initial data for 
the Einstein vacuum equations which contain trapped surfaces. The significance of this problem 
stems from the singularity theorems of General Relativity - or rather, the issues they raise. Take 
the Penrose theorem [1] for example. This shows that a spacetime N with reasonable matter 
(in particular: a vacuum spacetime) evolving from a non-compact (in particular: asymptotically 
flat) Cauchy surface M C N cannot be null-geodesically complete when M contains a trapped 
surface (TS). The theorem still holds when 'trapped surface' is replaced by 'outer trapped 
surface'. By a TS wc shall always mean an outer trapped surface. An obvious question to ask, 
is under what circumstances the TS-condition is actually satisfled. From the point of view of 
the initial-value problem this question really consists of two parts: 

1. Can one flnd conditions on M, e.g. in terms of 'free data', in order for M to contain a 
TS? 

2. What are conditions on M so that there is a TS in a Cauchy development of M? 

Thirdly, since the singular behaviour predicted by the singularity theorems is merely the exis- 
tence of incomplete causal geodesies, one would like to understand the local and global nature 
of these 'singularities'. In particular, one would like to know 

3. Does N satisfy some version of the cosmic censorship hypothesis (CCH)? 

Clearly, question 3.) is the most interesting one. Since 2.) and 3.) are by their nature ques- 
tions concerning the global Cauchy problem (see Chrusciel [2] for a review of what is known in 
this area) they are also very difficult. Only in the case of gravity coupled to a massless scalar 
field, and with spherical symmetry imposed on the geometry and the source, has a full analysis 
of these issues been achieved in the work of Christodoulou [3]. 

The present paper is devoted to the study, for initial data corresponding to pure gravitational 
radiation, question 1.) in the above list. Although this is a more modest task than attacking 
2.) or 3.), it is widely assumed that a full understanding of 1.) should also shed light on 3.). 
The reason for this lies in the belief, which is of course just one aspect of CCH, that the 
maximal Cauchy development of an asymptotically flat slice M with trapped surfaces is a black 
hole spacetime. In particular, singularities should be inside the apparent horizon which, in 
turn, is surrounded by an event horizon. No counterexamples are known to this 'event horizon 
conjecture' (see Israel [4] for steps twards a resolution). Initial data like the ones constructed in 
this paper are thus considered to be relatively 'safe' from the viewpoint of CCH. At the other 
extreme there have been studies of families of initial data the limiting members of which are 
singular, but without TS's. Such families, found numerically by Nakamura et al. [4], mimick the 
formation of naked singularities which Shapiro and Teukolsky [5] obtain in their simulations of 
collapsing localized dust. Since these solutions have only axial symmetry, the singularities are 
thought to be a much more serious threat to CCH than previously known counterexamples all 
of which are spherical. Sequences analogous to the ones in [4], but for the vacuum case, which 
are in a sense complementary to the vacuum sequences considered in our previous work [6], 
have recently been found numerically by Abrahams et al. [7] (see also O'Murchadha [8]). These 
may be taken as possible evidence for the violation of CCH even in the vacuum case, but the 
issue is completely open. 
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There has been considerable interest recently in conditions on Cauchy data in order to 
have trapped surfaces. Only the spherically symmetric case is fully understood so far [9]. For 
generalizations of these results, mostly motivated by Thome's Hoop Conjecture, sec Malcc [10] 
and Flannagan [11]. For a completely different approach, based on a suggestion by Penrose, see 
Tod [12] and Barrabes et al. [13]. (Compare, also, the results of Schoen and Yau [14], together 
with O'Murchadha [15].) 

In our work [6] we have shown the existence of time symmetric, asymptotically flat, vacuum 
initial data with TS's, thus settling a question raised by Bartnik [16]. (For earlier numerical 
evidence see Eppley [17] and Miyama [18].) Our construction is based in an essential manner 
on the conformal approach (see Choquet-Bruhat and York [19]) to solving the constraints of 
General Relativity. In the time-symmetric case this consists of picking a 3-mctric g on M which 
tends to a fixed flat metric 6 near infinity and seeking as the physical metric a metric g' in 
the conformal class defined by g, also tending to S at infinity. This gives rise to an elliptic 
equation for the conformal factor <I>, the Lichnerowicz equation, together with the boundary 
condition that $ tends to 1 at infinity. Some condition on g has to be satisfied in order for 
there to be a solution of this problem. This condition has the form of an inequality on a global 
conformal invariant associated with g, the Yamabe number Y (Lee and Parker [20], see also 
O'Murchadha [21]), namely 1" > 0. We now choose a sequence of metrics called a critical 
sequence (CS), for which Y > 0, but tending as n ^ oo to a metric for which y = 0. It 
follows that the ADM masses m„ increase without bound as n — > oo. It furthermore can be 
shown that higher-order multipoles diverge no more rapidly than the mass for large n. Prom 
this, in turn, one can infer that TS's form. The present work extends the above construction 
to the maximal (non-time symmetric) case, provided the suitably rcscalcd extrinsic curvature 
stays bounded as n increases. Then the ADM mass is dominated by its time symmetric part, 
so that essentially the same conclusion obtains as before. In the context of the discussion at the 
beginning of this section some features of the above scenario arc worth pointing out. Firstly, 
for large n, CS's can not result from time evolution via the Einstein equations. To see this, 
recall that maximal, asymptotically fiat slices in a single spacetime asymptotically have to be 
boosts or translations of each other (see Bartnik et al. [22]). We now invoke the fact that the 
spatial metrics g in our construction arc, in addition to being asymptotically flat at infinity to 
leading order, asymptotically spherical ( "Schwarzschildian" ) to the next order. This implies that 
asymptotic boosts are no longer allowed, and we are left with the possibility of slices which are 
mapped into each other by asymptotic translations. But then the mass m has to be conserved 
which contradicts limn^^ijin = oo. This suggests that one way to improve our results would 
be to try to establish a universal bound on Y, below which g' has TS's so that it would no 
longer be necessary to consider sequences with y — > 0. Secondly, it would be interesting to have 
a better understanding of Y in terms of conformal geometry on the local or quasilocal level and 
use this to study possible analogues of the Hoop Conjecture in the vacuum case. 

The critical sequences discussed in [6] are explicitly time-symmetric. This means that the 
TS's which form are both future and past trapped. This, in turn, means that the spacetimes 
evolving from these data for large n are both future- and past-incomplete. Thus, if cosmic 
censorship holds, we should have a configuration which starts out as a white hole in the distant 
past, expands to reach an instant of time-symmetry and then collapses to form a black hole in 
the distant future. In contrast, in our current work we break the time-symmetry assumption on 
the initial hypersurface by introducing some nonzero (albeit relatively small) extrinsic curvature. 
Thus, as we move along a critical sequence, we do not expect future and past apparent horizons 
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to appear simultaneously. Of course, if we move sufficiently far along any critical sequence the 
effects of the ever-increasing ADM mass will swamp the extrinsic curvature and we return to 
the time-symmetric situation with past as well as future singularities. Still, along any critical 
sequence we should have a band of data, solutions to the constraints, with a future apparent 
horizon but no past apparent horizon (if they appear in the opposite order, just multiply the 
extrinsic curvature by —1). These data now form possible candidates to generate spacetimes 
which represent the inflow of gravitational radiation from past null infinity of such an amount 
as to cause gravitational collapse and the formation of a black hole. We realize that the absence 
of a past horizon in an initial data set is no guarantee of regularity into the past (see question 
2.) of our first paragraph). However, it certainly is a necessary condition. Thus we can only 
claim that it is plausible that at least some of these initial data sets that are singular to the 
future, are regular to the past. 

Recall that in the conformal approach to solving the constraints [19] in the maximal case we 
have, in addition to the 'background' metric cjab^ ^ background extrinsic curvature Pah related 
to p^^, the physical extrinsic curvature, by = ^^'^^ah- In the calculation presented here 
we move the background geometry along a critical sequence while keeping the background 
extrinsic curvature bounded. The obvious alternative approach is to hold 'g fixed and scale 
up the extrinsic curvature p. While there are obvious technical difficulties we expect that this 
approach would also work and give initial data which are significantly more time-asymmetric 
than those produced here. This approach deserves further study. 

We restrict our attention to maximal initial data. The maximal condition allows a separa- 
tion of the momentum constraint from the Hamiltonian constraint and reduces the problem of 
constructing initial data from a quasi-linear coupled system to the problem of solving a single 
quasi-linear elliptic equation. We probably can handle the situation where the trace of the ex- 
trinsic curvature, while not vanishing, is small. Unfortunately, a systematic technique along the 
lines of this paper giving TS's in the non-maximal case requires a better understanding of the 
coupling between the Hamiltonian and momentum constraint than is presently available. 

Our paper is organized as follows. In Sect. 2 we recall the conformal approach to solving the 
Hamiltonian constraint in the maximal case and explain our basic assumptions. As in [6] we 
use a technique based on conformal compactification of 3-space M. On the resulting compact 
manifold M with metric g we use for simplicity, instead of the conformal invariant Y the quantity 
Ai, the lowest eigenvalue of the conformal Laplacian on {M,g). The latter is conformally non- 
invariant, but is positive (zero) iff Y is positive (zero), and this sign condition is the only one we 
use. In Sect. 3 the notion of critical sequence is introduced and the theorem on TS-formation 
is proved. In the Appendix we prove an existence and uniqueness theorem for the Lichnerowicz 
equation appropriate for our setting. 

2 Solving the Hamiltonian Constraint 

The basic formula underlying the standard [19] method of solving the Hamiltonian constraint 
is as follows. Define the conformal Laplacian associated with the 3-metric gab by 



Lg --Ag + lnlg], 



(2.1) 
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:= g^-^DaDt, being the standard Laplacian acting on scalars and TZlg] the scalar curvature 
of g. Then, if we put 

9ab = ^^9ab, $' = a;-i$ a; > 0, (2.2) 

we find that 

Lg><^' = (jj-^Lg^. (2.3) 
If in particular $ > and we set a; = we have = 1 and there results the identity 

^n[g'] = ^-'Lg<S>. (2.4) 

Thus TZ[g'] is zero if and only if <I> satisfies the 'linear Lichnerowicz equation' Lg^ = 0. Recall 
that 7^ = is just the initial-value constraint in the time-symmetric case. 

More generally, let us be given a trace-free, divergence-free w.r.t. to gab (TT-) tensor Pab, 

i.e. 

5«V = 0, D'^pab = 0. (2.5) 
Notice that Equ.-s (2.5) are conformally covariant if Pab is transformed like 

Pab = ^'"^Pab- (2.6) 
Suppose now that a positive scalar is subject to the equation 

Lg^ = * > 0, (2.7) 

then the pair (5^^ = ^^gab, p'ab = ^~^Pab) satisfies 

n[g'] = $-5$- W = 9"'''9"'PacP'bd (2.8) 

and 

D>'ab = 

that is to say the initial-value constraints of General Relativity in the maximal, source-free case. 

The set of pairs {gab 1 Pab) 1 where pab is TT, subject to the equivalence relation gab ^ g'ab 
Pab ~ p'ab '^ith g'^^^ = ijJ^gab and = uj^'^Pab for some w > 0, are called the 'free data' for 
vacuum G.R. The word 'free' here is misleading for two reasons. Firstly, Pab being TT comes 
from solving a system of partial differential equations. We will comment on how this can be 
done. Secondly, gab has to be such that the Lichnerowicz equation (2.7) has a positive solution 
<5 with the required global properties given the appropriate global properties of {M, gab, Pab)- 
To these we now turn. 

We want the 'physical' data 5^^, to be in a suitable sense asymptotically flat on some 
physical manifold M. Our idea here is to obtain {M,g',p') by 'conformal decompactification' 
from {M,g,p) where M is a compact manifold. Recall that it is often convenient to describe 
asymptotic fiatness of a Riemannian manifold (M,g') by the existence of a manifold M = 
M U {A}, A being a point, and a function on M, going to zero at A in such a manner that 
gab = ^"^g'ab extends in a sufficiently regular way to a metric on M. But this transition from M 
to M has a drawback: whereas in the physical picture we have metrics, in the unphysical picture, 
due to the gauge freedom in we only have conformal metrics, and these are more complicated 
to deal with. In the present situation, however, it is precisely the conformal structure associated 
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with gab which is among the basic objects, so the above freedom is in fact desirable. Thus we 
take M to be compact and gab a smooth metric thereon. Pick a point A G M, the point-at- 
infinity. (The results of the present section could easily be generalized to the case where there 
is an arbitrary but finite number of points-at-infinity.) We next define what we mean by an 
asymptotic distance function (ADF) w.r. to A. O^/^ is called an ADF, provided CI is smooth on 
M, positive outside A and satisfies (Oa := Da^, 0,ab '■= DaDb^, a.s.o.) 

^1^ = 0, ^^a|A = 0, {nab-2gab)\jy = 0, f^afccU = 0. (2.9) 

For example, could be taken to be the square of the geodesic distance from A. Given gab 
and O, the metric gab, defined by gab = ^~^9ab is easily seen to be asymptotically flat on M 
and have vanishing ADM mass. More precisely, if are coordinates centered at A, in which 
9ab\A ~ ^ab and dgablj^ = 0, then = Q~^x°' are coordinates in a punctured neighbourhood of 
A, in which 

gab = Sab + 0{n), (2.10) 

dab being theKronecker delta. Thus in order for ^'^gab to be an asymptotically fiat metric the 
correct boundary condition for $ is something like 

$ = J^-V2 + ^ + 0{n^/^), m = constant. (2.11') 
In fact, we shall have to weaken (2.11') by only requiring that 

$ = J^-V2 + ^ + 0(J^£/2-)^ 0<£<1, (2.11) 

9$ = 9(0-1/2) ^0(^(£-i)/2)_ (2.12) 

When $ obeys (2.11,12), then ^f^;, = $^5a6) in the x"-coordinates, satisfies the asymptotic 
relations 



9ab 



l + ^Ua6 + 0(0(^-+i)/2), (2.13) 



dg'ab = dil + ^— j 6ab + 0{n^/^). (2.14) 
Writing = SabX^-x^, we have 

Q = + 0{r~^) near A (2.15) 

so (2.13,14) are in fact the standard conditions for an asymptotically Schwarzschildian 3-metric. 
The boundary conditions for p'^^^, formulated in terms of Pa6) are as follows: 

Pab = 0{n-'/^+'/% (2.16) 
dpab = 0(J^-lV4+./4) 

which implies that p^;, = ^~^Pab satisfies 

p'^^ = 0(0^/4+-/^) = 0(r-=^/2-^/2), (2.17) 
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The connection of the above 'unphysical' formulation with thejtandard physical formulation 
of the Lichnerowicz equation is as follows. We define fields on M 

$ = = (2.18) 

so that $ — > 1 at infinity and Pab has the same asymptotic behaviour as (Equ. (2.17)). 
Furthermore 

^$ = l^-^p^^jf^b ^. (2.19) 

The condition that the extrinsic curvature falls off faster than r^-^/^ ig exactly the one that 
guarantees that the ADM mass be finite and well-defined [24]. Naively, to leading order the 
physical Lichnerowicz equation looks like A$ + |pa6P"^ = with PabP""^ = 0{r~^~^). Hence 
$ = 1 +m/2r + 0(r-i-^). This gives us the standard Schwarzschildian behaviour at infinity. 

Using the results of Chaljub-Simon [23] one can construct many PabS satisfying (2.5) and 
the boundary conditions (2.16) in the case where M = (M = R^). 

We now come to the issue of existence of solutions $ > 0, satisfying 

Lg^ = l^-^Pabp''^ onM\A (2.20) 

o 

together with (2.11,12). 

First recall that Lg is an essentially self-adjoint operator on C°°(M) w.r. to the standard 
L^-inner product associated with the metric g. In addition, the spectrum of Lg is a pure point 
spectrum and is bounded from below. Let Xi{g) be its lowest eigenvalue. 

Theorem: Let Xi{g) > 0. Then (2.20), (2.11,12) have a unique solution which is smooth on 
M\A. 

Proof: It is known (see e.g. [20]) that the problem can be uniquely solved when pab is zero. 
Call the corresponding solution G and write 

^ = G + h. (2.21) 

Then (2.20) takes the form 

Lgh = ^(1 + G-'h)-'G-'pabP''', (2.22) 

where G~^h simply means pointwise multiplication. Note that the r.h. side of (2.22) diverges 
at the point A, but does so in a relatively mild fashion, i.e. like 0{r~'^~^^), < £ < 1. This is 
enough for the existence theorem proven in the Appendix to hold. 

3 Critical Sequences 

Similarly as in [6] we call a sequence gn of 3-metrics on M with Xi{gn) > critical when it 
converges uniformly together with all its derivatives to a metric g^o with Xi{goo) = 0. We also 
require a critical sequence (CS) to be such that the metrics gn and the connections dg-n all 
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coincide at the point A, so that an ADF fi^/^ can be found which is independent of n. We 
furthermore assume that for each n we are given a TT-tcnsor pab which is smooth in M \ A and 
satisfies (2.16), where the constants involved in the O-symbols are understood to be independent 
of n. The existence of such a sequence of Pafe's, in the case where M = S^, follows from the work 
of Chaljub-Simon [23]. (Basically the reason for this is that the criticality of the limiting metric 
^oo plays no role for the existence of a TT-tensor relative to goo-) For typographical reason we 
omit the index "n" in the following expressions. Since Xi{g) > 0, the operator Lg has a Green 
'function' G{x,x'). The function G{x) solving the time-symmetric constraints is, of course, 

G{x) = G{A,x). (3.1) 

Thus the solution to Equ. (2.22) can be written as 

h{x) = \[ G{x,x')[l + {G-^h){x')]-\G-'pabP''''){x')dVg{x'). (3.2) 
o Jm 

Let Q}/'^{x^ x') be the generalization of an ADF, with the point A replaced by an arbitrary point 
x' G M. (Note that f2(a;,x') depends on n.) Then 

G(x, x') = n-'/^{x, x') + ^ + 0(J^V2). (3.3) 

But 

n{x, x') = n{x', x) [1 + o{n)] (3.4) 

and, by virtue of the self-adjointness of Lg, 

G{x,x') = G{x',x), (3.5) 

so that n{x') in (3.3) is in fact constant and thus identical with the ADM mass of the time- 
symmetric contribution to the initial-data set. Next, calling 

G{x, x') - n-^/^{x, x') = ^{x, x'), (3.6) 

it follows from the results in [6], that lim„_+ooA*n = oo and there are positive, n-independent 
constants C, C such that 

C/i„ <ij< C'fln. (3.7) 

(This is proved in [6] for the Green function G{x), i.e. for x' at A where g, dg are kept fixed. It 
is not difficult to see that the statement remains true in the present situation.) 
Furthermore, in an arbitrary chart centered at x, 

D\x- x'\'^ <n{x,x') <D'\x- x'l"^ (3.8) 

where I ^ and £) > 0, D' > 0, both independent of n. Thus 

1 



G{x,x') < const 

and 



X — X' 



j7 + IJ'n const 



(3.9) 



const 



1 

■j— ^ -I- const IJ,n 



< G{x'). (3.10) 
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Similarly, 



dG{x, x') < const 



,T — X 



'12 



+ jJLn const 



(3.11) 



In order to estimate the integral in Equ. (3.2) for large n, it suffices to take the 'worst' part, 
i.e. the contribution where the integration is carried out over a neighbourhood A^a of A and 
where x lies in that neighbourhood (so that both the singularity of G(x, x') aX, x = x' and the 

singularity of {G~'^ PahP'^^){^') at x' = A contribute). Inserting (3.9,10,11) together with (2.16) 
into (3.2) and using that the volume element dVg{x') has a uniform bound in terms of the 
Euclidean volume element, we are left with an integral of the form 



= [ 

Jr'< 



r'<R \ \X — X' 



j:+fin const) + Hn const) r'~^'^^r''^dr' sin e'de' dip', (3.12) 



r < R, < £ < 1, 



where r' , 6', ip' are to be thought of as polar coordinates coming from a chart x' centered at A. 
After elementary integrations we find that 



0(/x- 



< e < 1. 



(3.13) 



Similarly, to estimate dh, it suffices, using (3.9,10,11) and (2.16), to consider the quantity 



Jn{x) 



1 



lr'<R l\X - X 

We easily see that 



'12 



+ jin const 



(1 + Hn const r') 



j-^dr' sin e'de' dp', Q < r < R. (3.14) 



Jjx) = Oir 



< £ < 1 



Recalling that $(a;) = G{x) + h{x) and using (3.7, 3.13, 3.15), it follows that 

$ = 0-V2 + ^+^„0(j^i/2) + o(^-.) 



(3.15) 

(3.16), 
(3.17) 



= d{^-^/^) + /i„o(i) + o(^7(^-i)/2)^ 

where $7^/^ is an ADF corresponding to the point A. Now consider level surfaces of il. When VL 
is sufficiently small but positive, these are embedded two-spheres. Let n'" be the unit normal 
with respect to the physical metric g'^^^ = ^^gab, pointing towards A. Then the surface is a TS, 
provided the quantity H defined by 



H = D',n"'+p',,n"^n"' 
is negative on the surface. We rewrite (3.28) in terms of unbarred quantities. Using 



we obtain for the first term in (3.18) the expression 



(3.18) 



(3.19) 



' „'a 



(3.20) 
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Prom the definition (2.9) of ft and the properties of a CS, we easily find 



- TTT^ ]nab = 4: + 0{n) uniformly in n. 



(3.21) 



Using (3.16,17), there results 



D'an 



3 r 



(3.22) 



We now have to estimate the p^^-term in Equ. (3.18). Using p'^f^ = and (3.19) we obtain 



1 $-3^^1/2 



ar>b 



(3.23) 



Due to (3.16) we have that <1> > O"^/^ for large n. Using this fact and (2.16) we see that 
the bracket in (3.23) is 0(1^^^"'"'^)/^). Combining this estimate with (3.22) it follows that H is 
negative provided the quantity 



< £ < 1 



(3.24) 



is negative, where A, B are suitable positive constants independent of n. It can be shown that 
there exist positive constants (independent of n) a and P such that 77 is negative for 



(3.25) 



when n is large enough. This completes the proof of our theorem. 



Appendix 

We do not attempt to formulate a theorem under the most general conditions. The equation 
we want to solve is of the following form 

Lgh={l + G-^h)-^p on M' = M\ A (A.l) 

where is a positive C'^-function on M. The source p{x) is also positive on M' and smooth. 
In coordinates x°' centered at A, it behaves like 

p{x) = 0(r-2+^), < £ < 1, (A.2) 

dp{x) = 0(r-^+^). (A.3) 

The metric g on M is such that \i{g) > 0. Hence the operator Lg in (A.l) has an everywhere 
positive Green function G{x,x'). Instead of (A.l), we first try to solve equation 

h = J^[h], (A.4) 
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where J^[h] is defined by 

{J^[h]){x)= I G{x,x'){l + G-^h)-\x')p{x')dVg{x'). (A.5) 
Jm 

We are thus seeking a fixed point of the mapping J^. When x' is in a local neighbourhood of x, 
G{x,x') satisfies 

G(x.x') = j^^-{-0(l), (A.6) 

with \x — x'p = gab{x){x — x')"'{x — x')^ . (A.6) and the fact that G > imply that T maps 
continuous, non-negative functions into themselves. Let -0 be the continuous function on M 
defined by 

= J^[0] (A.8) 

and U the space 

U = {h^C^{M)\{i<h{x)<^{x)} (A.9) 

which is clearly a closed, convex subset of the Banach space C^{M). Since J- is pointwise 
decreasing, it maps U into itself. Let h be in the image of U under Then, by (A. 7), it is 
easily seen to be in (M \ A) with a singularity at A of the form 

dh = 0{r'-^). (A.IO) 

From this one easily deduces that h is bounded in C^'^~^{M). Thus J'{U) is pre-compact in 
C^{M). Furthermore, for /ii,/i2 € U, there holds 

I^N-^NI <c|^i-^2|, (A.ii) 

where C is independent of hi, /i2- Thus !F, defined on J7, is continuous. By the Schauder fixed 
point theorem (see e.g. Gilbarg-Trudinger [25], p. 221) it now follows that there exists h £ U, 
such that 

h = T[h]. (A.12) 

Again, this implies that h G C'^'^^^(Af). Thus the function (1 + G~^h)~'^ in (A.5) is Lipschitz 
continuous. Furthermore, by virtue of (A. 2, 3), p satisfies 

\p{x) - p{y)\ < . jfr,^L, (A.13) 
IPV J Hyyji - [inin(|a;|, |2/|)]3-= ^ ^ 

near A. Using these facts in (A.5) together with standard potential theoretic estimates analogous 
to [25, p. 54] we infer that h G C^{M \ A) and satisfies (A.l). Furthermore ^ = G + h obeys 
the boundary conditions (2.11,12). 

Finally, for completeness sake, we show, along the lines of Cantor [26], that the above 
solution is unique. Suppose, assuming the contrary, that there are two solutions, $i and $2, 
solving Equ. (2.20). Then, defining 

9ab = G^gab, ^i = G-'^i, i = l,2, (A.14) 
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Lg$ = -Ag$ = G-V(ar, on M (A.15) 

where cr is a nowhere increasing function of #. Thus, with -0 defined by V' = — ^2, we have 

-AV' = o-(x,$i)-(7(x,$2)- (A.16) 

Suppose t/j is anywhere non-zero, say positive. Then there is an open set O = {x E M\il) > 0} 
which is not the zero-set with ^ = on dO. Thus 

AV' > (A.17) 

which, by virtue of the maximum principle, is a contradiction. Exchanging the role of $1 and 
$2 we find that, in fact, $1 = $2- This ends the uniqueness proof. 
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